Wavelet analysis of event by event fluctuations 



O 

o 



p. Manimarant and Prasanta K. Panigrahi^ ^ 

School of Physics, University of Hyderabad, Hyderabad, 500 046, India 
* Physical Research Laboratory, Ahmedabad, 380 009, India 
(Dated: February 2, 2008) 

The temporal fluctuations of produced hadron density in heavy ion collisions, modelled by 2D 
Ising model at temperatures Tc and below, are studied through a recently developed wavelet based 
fluctuation analysis method. At Tc, long-range correlated multifractal behavior, matching with the 
recently observed Hurst exponent iif ~ 1, is found. Below Tc uncorrelated monofractal behavior 
is seen. The correlation behavior compares well with the results obtained from continuous wavelet 
based average wavelet co-efflcient method, as well as with Fourier power spectral analysis. 
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1. INTRODUCTION 

The study of correlations and scaling behavior is an 
area of active research IJ. Various methods of anal- 
ysis like rescaled range analysis 2], structure func- 
tion 0, wavelet transform modulus maxima 0], de- 
trended fluctuation analysis and its variants are used to 
study the correlation behavior and fractal characteristics 
[ililllllillillllil- Recently, we have developed 
a new method based on discrete wavelet transform to 
study the scaling properties of non-stationary time series 
[T3l |. This procedure economically extracts the multi- 
fractal behavior in a time series. It has been applied to 
characterize current and voltage fluctuations in tokamak 
plasma and financial time series p^ . 

The goal of the present work is to apply this discrete 
wavelet based approach for the study of event by event 
fluctuations of hadron density produced during phase 
transition. The fact that, our method is a local one, 
makes it ideal for the analysis of fluctuations in a non- 
stationary data. The primary motivation for this work 
is the recent study of the event by event fluctuations of 
hadron multiplicities, carried out by Qin and Ta-chung, 
[isl IT^ using rescaled range analysis, which has shown 
scaling behavior. The Hurst exponent H , has been found 
to be 1. The geometry of the hadronization scenario in 
the context of heavy ion collisions, with two flavors, has 
led to the description of the same through a model like 
two dimensional (2D) Ising model. Hwa and co-workers 
have analyzed the scaling behavior of hadronic fluctua- 
tions through various methods like moment analysis etc., 
in the context of the above model |r7l lllll9ll2a | . 

A large number of studies have been done to analyze 
the characteristic behavior of hadron formation during 
quark hadron phase transition. Earlier studies regard- 
ing the fluctuations in particle production, like corre- 
lation analysis of hadronic fluctuations have been car- 
ried out using cluster expansion technique of Ursell and 
Mayer [23|. Bialas and Peschanski have studied multi- 
particle production during quark hadron phase transition 
through factorial moment analysis jlj . Using statistical 
methods, Ludlam and Slansky studied the mechanism of 
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FIG. 1: The integrated time series of average hadron densities 
after subtraction of the mean, [a] at T = 1.0, below Tc and 
[b] Tc — 2.27. The difference in behavior of the fluctuations 
at different temperatures is clearly visible. 



jarticle production during quark hadron phase transition 



Following Hwa et. al. jlTj, we have carried out 
the Monte-Carlo simulation of 2D Ising model to mimic 
the multi-particle production during quark hadron phase 
transition. The time series of event by event fluctua- 
tions of average hadron densities below and at phase 
transition have been measured for characterization. The 
spatio-temporal fluctuations of produced hadrons show 
correlation behavior which have been characterized us- 
ing wavelet transform techniques. We have made use 
of discrete and continuous wavelet transforms to analyze 
the hadronic fluctuations at various temperatures below 
and at critical temperature. We also study the scaling 
behavior through power spectral analysis. 

Wavelet transform is a mathematical tool which finds 

flications in various fields, ranging from finance 
biology [2^ [2^, to chemistry (2^, and physics 
^ [23 etc. Wavelet transform decomposes any given 
function or data in a manner, which achieves optimal 
time-frequency localization. Wavelets are classified into 
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discrete and continuous ones [30|, |3l|. We have made 
use of both the methods to analyze the fluctuations in 
this paper. In case of discrete wavelet analysis, we have 
used an approach which has been developed recently |l3j| . 
Average wavelet coefficient method, a continuous wavelet 
approach as well as Fourier analysis are also used to study 
the correlation behavior of the fluctuations and corrobo- 
rate our findings. 

This paper is organized as follows. Section II deals 
with hadron density as computed from the Ising model. 
In Section III, discrete wavelet analysis of time series 
has been carried out for extracting multifractal behavior. 
In section IV, continuous wavelet analysis and Fourier 
power spectral analysis have been implemented to sup- 
port the previous discrete wavelet transform approach. 
Section V concludes with discussions and future direc- 
tions of work. 



2. HADRON DENSITY AND THE ISING 
MODEL 

In the context of heavy ion collisions, a two dimen- 
sional geometry can be invoked, where 2D Ising model 
captures the dynamics of the phase transition in a two 
flavor scenario. Earlier studies by Hwa et.al., 0| de- 
scribe in detail the relationship between hadron density 
and the 2D Ising model. In this paper, we adopt the same 
procedure and deal with the correlation analysis and frac- 
tal characteristics through wavelet transforms. We have 
studied the hadronic fluctuations using both discrete and 
continuous wavelet transform as a time series analysis. 

The two dimensional Ising model is deflned by the 
Hamiltonian, 



(1) 



where s takes the values ±1 and the sum is over the 
nearest neighbors. In our Monte-Carlo simulation of 2D 
Ising model ^3^, we use Wolff single cluster algorithm 
[33] for spin flipping in a conflguration of a lattice size 
L = 256. We define the cell size to be with e = 4 and 
the whole lattice is subdivided into I = {L/e)'^ cells. The 
hadron density at a particular cell is defined by 



(2) 



where Ci is the net spin at cell i defined to be positive 
along the overall magnetization, i.e.. 



sign 




(3) 



Ai is the cell block of size located at i, 8 is the Heavi- 
side function and A is a constant. During the hadroniza- 
tion process q fluctuates from cell to cell, which reflects 



in the average hadron density in event to event. The 
average density is calculated as; 



e=l i=l 



(4) 



Here, i denotes the cells and e is the configurations cap- 
tured, where iV is the number of events simulated. We 
have calculated the equilibrium time series of the average 
hadron density over iV = 5 x 10^ iterations. 



3. WAVELET BASED FLUCTUATION 
ANALYSIS 

Considering the time series of average hadron density 
{{p)i) from i = 1, N , we capture configurations at var- 
ious temperatures like below, and at critical temperature. 
In most applications, the index i will correspond to the 
time of the measurements. The profile (1^) of the time 
series is obtained by subtracting the mean and taking 
cumulative sum of the average hadron density. The pro- 
files have been analyzed by the wavelet based fluctuation 
analysis method to study the scaling behavior (l^j. We 
calculate the scaling exponents h{q) for various moments 
q, in this case q varies from —10 to 10. For q = 2, h{q) is 
the Hurst exponent, which is one of the measures of cor- 
relation behavior in the time series analysis: < < 1. 
For persistent time scries H > 0.5, H ~ 0.5 for uncor- 
related time series and H < 0.5 for anti-persistent be- 
havior. Recently, the anti-persistent regime < H < 0.5 
has been mapped to A3 statistics, widely employed in 
random matrix theory 34]. The main aim here is to 
study the fractal characteristics of the time series and its 
multifractal nature. For the time series possessing multi- 
fractal behavior, h{q) decreases with increasing values of 
q. For mono-fractal time series h{q) is constant for all 
values of q. Fig. 1 depicts the proflle of the time series 
of event by event fluctuations of average hadron density, 
[a] below T^, and [b] at T^. 

We use a discrete wavelet based fluctuation analysis 
method (which is analogous to multifractal detrended 
fluctuation analysis method) to study the existence of 
scaling behavior of the time series. Our earlier paper 
gives the detail procedure of this approach The 
fluctuations are extracted using a wavelet (Db-8) from 
Daubechies family. The fluctuation function Fq{s) is ob- 
tained for various scales s. The power law manifests itself 
as a straight line in the log-log plot of Fq{s) versus s, for 
each value of q. 



(5) 



From the analysis of hadron density time series, below Tc, 
uncorrelated monofractal behavior and long-range corre- 
lated multifractal behavior at Tc are found. These results 
are shown in Fig. 2. It was found that the Hurst expo- 
nent H — h{q = 2) ~ 1 matches with that of Qin et. al. 
[? ]. Very interestingly, multifractal behavior is clearly 
seen for event by event fluctuations at Tc- 
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FIG. 2: For time series [a] Below Tc, h{q) values shows lin- 
ear behavior for different values q indicating the monofractal 
behavior and [b] the non-linear behavior of h{q) values for 
different values of g at Tc, shows clearly the long-range corre- 
lation and multifractal nature. 



4. AVERAGE WAVELET COEFFICIENT 
METHOD 

We now make use of average wavelet coefficient mettiod 
to corroborate the above findings. This is a continuous 
wavelet based approach, which has been used to study 
the temporal correlations of the fluctuations in time se- 
ries. We obtain the Hurst measure by transforming the 
time series, f{t) into wavelet domain. The continuous 
wavelet transform of a function f{x) is given by 



FIG. 3: Log-Log plot of average wavelet coefficients VK[/](a) 
versus scale a, indicates [a] below Tc, uncorrelated behavior 
{H ~ 0.5) and [b] at Tc, long-range correlation {H ~ 1). 



o 

- 


slope a = 2 


' [a] 





where the mother wavelet ■0a, 6 is defined as, 



V'a,6(a;) = 



a 



(6) 



(7) 



Here a is the scaling parameter (a > 0) and 6 is the 
translation parameter (— oo < b < oo) and ip*{x) is the 
complex conjugate of i^{x). 

It has been shown that ISal , 



W[n{a) = {\W[f]{aM)^^-^ 



+H 



(8) 



At a given scale, the wavelet energy of all locations 
have been averaged yielding M^[/](a). The same averag- 
ing is followed for all scales. Thus for a self-affine time 
series, the scaling exponent (| -I- H) is measured from 
the slope of the log-log plot of T4^[/](a) versus scale a. 
The results are shown in Fig. 3. They clearly indicate 
that the fluctuations possess self-affine nature and the 
Hurst exponent is the measure of correlation. From the 
obtained H , we found that the hadronic fluctuations are 
uncorrelated below Tc {H = 0.5) and at Tc it possesses 



FIG. 4: Log-Log plot of Fourier power spectrum P{s) versus 
frequency s, shows [a] below Tc, Brownian motion behavior 
a = 2 and [bl at Tc, fractional Brownian motion a = 3. 



long-range correlation (H ~ 1). This matches with our 
previous discrete wavelet analysis. 

We have also analyzed the scaling behavior through 
Fourier power spectral analysis. 



Yit) exp{-2mst) dt 



(9) 



Here Y(t) is the accumulated fluctuations after subtract- 
ing the mean {¥). It is well known that, P{s) ^ s^". 
The obtained scaling exponent a is compared with Hurst 
exponent by the relation a — 2H + 1. For the time series 
far below from Tc, the scaling exponent a = 2 which re- 
veals uncorrelated Brownian motion behavior and at Tc, 
the scaling exponent a = 3, which reveals long range cor- 
related fractional Brownian motion behavior as shown in 
Fig. 4. 
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5. CONCLUSION 

We have found that the event by event fluctuations of 
hadron density, as modelled by 2D Ising spins, possess 
long-range correlated multifractal and fractional Brown- 
ian motion behavior at Tc- This is due to the correlation 
arising from phase transition. Uncorrelated monofractal 
and Brownian motion behavior below Tc is observed as 



per physical expectations. The analysis has been carried 
out through discrete wavelet wavelet based fluctuation 
analysis method, which is well suited for non-stationary 
data. The continuous wavelet based average wavelet co- 
efficient method and Fourier analysis fully supported the 
above conclusion. 

We thank Prof. J. C. Parikh and M. S. Santhanam for 
useful discussions. 
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